Nonequilibrium dynamics of a spin-3/2 Blume Capel model with quenched random 

crystal field 
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The relaxation and complex magnetic susceptibility treatments of a spin-3/2 Blume-Capel model 
with quenched random crystal field on a two dimensional square lattice are investigated by a method 
combining the statistical equilibrium theory and the thermodynamics of linear irreversible processes. 
Generalized force and flux are defined in irreversible thermodynamics limit. The kinetic equation 
for the magnetization is obtained by using linear response theory. Temperature and also crystal field 
dependencies of the relaxation time are obtained in the vicinity of phase transition points. We found 
that the relaxation time exhibits divergent treatment near the order-disorder phase transition point 
as well as near the isolated critical point whereas it displays cusp behavior near the first order phase 
transition point. In addition, much effort has been devoted to investigation of complex magnetic 
susceptibility response of the system to changing applied field frequencies and it is observed that 
the considered disordered magnetic system exhibits unusual and interesting behaviors. Furthermore, 
dynamical mean field critical exponents for the relaxation time and complex magnetic susceptibility 
are calculated in order to formulate the critical behavior of the system. Finally, a comparison of 
our observations with those of recently published studies is represented and it is shown that there 
exists a qualitatively good agreement. 

PACS numbers: 64.60.Ht, 75.10.Hk, 75.30.Cr, 76.60.-k 

Keywords: Nonequilibrium phase transitions, Onsager theory of irreversible thermodynamics, Quenched 
random crystal field. 



I. INTRODUCTION 

Investigation of disorder effects on the critical phe- 
nomena has a long history and there have been a great 
many of theoretical studies focused on disordered mag- 
netic materials with quenched randomness where the ran- 
dom variables of a magnetic system such as random fields 
[l|, [|| or random bonds [3, Hj may not change its value 
over time. On the other hand, quenched crystal field di- 
luted ferromagnets constitute another example of mag- 
netic systems and equilibrium properties of this type 
models have been investigated under multifarious ap- 
proximation methods [^-Q. Besides, effects of disorder 
on magnetic systems have been systematically studied, 
not only for theoretical interests but also for the identifi- 
cations with experimental realizations ■ It has been 
shown by renormalization group arguments that first- 
order transitions are replaced by continuous transitions, 
consequently tricritical points and critical end points are 
depressed in temperature, and a finite amount of disorder 
will suppress them 

Since the time in which Ising model j_13i] was invented 
there exists a limited number of studies including the dy- 
namic nature of the system under the influence a small 
perturbation. Even though the physical investigations 
regarding these systems bring about a lot of mathemat- 
ical difficulties, the nonequilibrium systems are in the 
focus of scientists because they have an unusual and in- 
teresting dynamic behavior. It is known that examina- 



tion of physical properties of the nonequilibrium systems 
allows us to analyze various dynamical concepts. One of 
them is the relaxation time of a considered system which 
exhibits a divergence near the critical and multicritical 
points and gives remarkable information about dynamic 
properties. As far as we know, on the theoretical pic- 
ture, the relaxation behavior of a limited number of sys- 
tems has been investigated by making use of Onsager 
reciprocity theorem [lj, EH ■ This t ype of investigation 
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was first performed by Tanaka et al. [16l | near the contin- 
uous phase transition point for the spin- 1/2 Ising model 
by means of Onsager reciprocity theorem. Next, AB- 
type ferromagnetic and antiferromagnetic models have 
been investigated by Barry [17| . Barry and Harrington 
fl~8| using the similar formulation, respectively. Addition- 
ally, in Ref. [Tc| . dynamic initial parallel susceptibility 
expression is analyzed near the Neel temperature for con- 
sidered applied field frequencies. Recently, spin-1 Ising 
model has been analyzed by Erdem and Keskin within 
the framework of the entropy production [l|| [2(| and it 
is found that one of the relaxation times diverges at the 
continuous phase transition point while the other relax- 
ation time remains finite. The relaxation process of the 
pure spin-3/2 Blume- Emery-Griffiths (BEG) model with 
bilinear and biquadratic exchange interactions has been 
analyzed near the critical 2l| and multicritical points 
22]. In addition, thermal variations of the relaxation 
time of a metamagnctic Ising model has been investi- 
gated by Gulpinar et al. [23f . Very recently, following 
the same methodology, spin-1 Blume-Capel (BC) model 
with quenched random crystal field has been studied in 
detail and it is observed that the relaxation time has 
a jump discontinuity at the first-order phase transition 
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point whereas it diverges in the neighborhood of multi- 
critical points such as tricritical, critical end point and 
bicritical end points [24]. Moreover, from the experimen- 
tal point of view, a great many of experimental studies 
have been dedicated to the better understanding of dy- 
namic critical phenomena in magnetic systems |25l - l29| . 

On the other side, it is known that the dynamic or 
ac susceptibility measurement, in which a time varying 
oscillating magnetic held is carried out to a sample, is a 
powerful method for analyzing the dynamic evolution of 
the considered real magnetic system. It is obtained from 
the dynamic response of the system to time dependent 
magnetic field and many studies have been performed 
regarding the magnetic relaxation of cooperatively inter- 
acting systems, such as nanoparticles |30l - l32| . spin glasses 
[33j , high T c systems [34[ and magnetic fluids [35| . From 
the theoretical point of view, in Ref. [36] , Acharyya and 
Chakrabarti probed thoroughly the real and imaginary 
parts of magnetic susceptibility in the neighborhood of 
phase transition temperature of a spin- 1/2 Ising model 
under the time dependent oscillating external magnetic 
field by using Monte Carlo (MC) simulation with peri- 
odic boundary conditions. Additionaly, Erdem investi- 
gated the magnetic relaxation in a spin-1 Ising model 
near the second-order phase transition point where time 
derivatives of the dipolar and quadrupolar order param- 
eters are treated as fluxes conjugate to their appropriate 
generalized forces in the sense of irreversible thermody- 
namics 1371 ] . next, the same author has analyzed the fre- 



quency dependence of the complex susceptibility for the 
same system [38|]. In addition to these, with the helping 
of Nelson's method, a systematic investigation contain- 
ing frequency, momentum and temperature dependent 
response function has been achieved for Ising-type system 
below the critical temperature [39j ] . Very recently, a com- 
prehensive study including the equilibrium and nonequi- 
librium antiferromagnetic and ferromagnetic susceptibil- 
ities of a metamagnetic Ising system in the vicinity of 
order-disorder transition point has been done by benefit- 
ing from mean field approximation (MFA) |40j . 

Apart from these, the effects of impurities on driving- 
rate-dependent energy loss in a ferromagnet under the 
time dependent magnetic field have been analyzed by 
Zheng and Li [41] by using several well defined mod- 
els within the frameworks of MFA and MC, and they 
found using MFA that, the hysteresis loop area is a power 
law function of the linear driving rate as A — Ao oc ft,* 3 , 
where, Ao,h and /3 are the static hysteresis loop area, 
the linear driving rate and scaling exponent of the sys- 
tem, respectively. Very recently, the quenched site and 
bond diluted kinetic Ising models under the influence of 
a time dependent oscillating magnetic field have been 
analyzed by making use of effective field theory 4^, 43[ 
on a two dimensional honeycomb lattice. In Ref. [42j, 
the global phase diagrams including the reentrant phase 
transitions are presented by the authors for site diluted 
kinetic Ising model, and they showed that the coexis- 
tence regions disappear for sufficiently weak dilution of 



lattice sites. Following the same methodology, the au- 
thors have concentrated on the influences of quenched 
bond dilution process on the dynamic behavior of the 
system. After some detailed analysis, it has been found 
that the impurities in a bond diluted kinetic Ising model 
give rise to a number of interesting and unusual phenom- 
ena such as reentrant phenomena and the impurities to 
have a tendency to destruct the first-order transitions 
and the dynamic tricritical point (43| . Furthermore, it 
has also been shown that dynamically ordered phase re- 
gions get expanded with decreasing amplitude which is 
more evident at low frequencies. 

As far as we know, the nonequilibrium properties of the 
quenched disordered systems such as crystal field (and 
also bond or site) diluted ferromagnets under the influ- 
ence of a small magnetic field perturbation have not yet 
been investigated for higher spin models. These types 
of disorder effects constitute an important role in the 
real magnetic material science, since the quenched disor- 
der effects may induce some important macroscopic in- 
fluences on the material. Therefore, we believe that the 
investigation of the effects of quenched randomness on 
the nonequilibrium properties of the Ising model and its 
derivations still need particular attention. Hence, in this 
work based on MFA, we intend to study the nonequi- 
librium properties of a spin-3/2 BC model by introduc- 
ing the quenched crystal field dilution effects. Here, we 
used a method combining statistical equilibrium theory 
and the thermodynamics of irreversible processes to ana- 
lyze the magnetic relaxation behavior and dynamic com- 
plex susceptibility of a spin-3/2 BC model with quenched 
crystal field randomness. In this context, by means of 
the MFA to obtain magnetic Gibbs free-energy, a gen- 
eralized force and a generalized current are determined 
within the framework of irreversible thermodynamics in 
the neighborhood of equilibrium states. After that, the 
kinetic equation for the time dependent magnetization 
is obtained and a relaxation time is derived and also 
thermal and crystal field variations are examined in the 
neighborhood of critical and isolated critical points. In 
addition, we calculated the complex magnetic suscepti- 
bility expression and analyzed near the isolated critical 
as well as second order phase transition points for se- 
lected Hamiltonian parameters. It can be said that the 
considered disordered magnetic system exhibits unusual 
and interesting behaviors. The aforementioned features 
will be discussed in detail in later sections. 

The remaining part of the paper is as follows: In Sec. 
im we briefly present the formulations, the results and 
discussions are presented in Sec. IIII1 and finally Sec. HVl 
contains our conclusions. 



II. FORMULATION 

In this section, we give the formulation of present study 
for a two-dimensional square lattice. The Hamiltonian of 
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the spin-3/2 BC model is given by: 



H 



(<i> 



(1) 



where the first term is a summation over the nearest 
neighbor spins with Si — ±3/2, ±1/2 and A<, in the 
second term, represents random crystal field which is 
distributed according to a given probability distribution 
function: 

G(Ai) = 1 (5(Ai - A(l + a)) + 5(Ai - A(l - a))) , 

(2) 

here a > 0, and half of the lattice sites subject to a crys- 
tal field A(l + a) and the remaining lattice sites have a 
crystal field A(l— a). In order to examine the nonequilib- 
rium properties of the considered magnetic system, first 
we have to focus on the equilibrium properties following 
the definition given in Ref. In this context, we have 
derived the Gibbs free energy (G = U — TS — BNm) by 
making use of MFA in the following form: 



g(m,t,d,b) = G/N 
-| [log (2exp £ 3 + 2 exp (-%) 4)] (3) 



[log (2 exp £ 3 + 2 exp (-^) 4)] 



where z, t = T/J, b = B/J and d = A/ J are the 
coordination number, the reduced temperature, the re- 
duced magnetic field and reduced crystal field, respec- 
tively. In addition, for the sake of simplicity, the coeffi- 
cients = 1, ..,4) are defined as follows: 



p _ d(l-a) 



cosh ( 



(4) 
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In the absence of the external applied field (b) and by 
benefiting from the minimization condition of the free 

energy at the equilibrium = ) , the mean field equa- 

tion of state is found as follows: 



34 + e 2£l 4 34 + e % 4 



44±4e 2£ i4 44 + 4e^4 



(5) 



here £ 5 = sinh(2|^) and £ 6 = sinh(^). This equa- 
tion may be solved without difficultly by an iterative 
procedure, i.e. Newton-Raphson method. Because the 
detailed equilibrium investigations of the considered sys- 
tem are analyzed in Ref. [44j], we shall only give a brief 
summary here as follows: Equilibrium phase diagrams of 
the spin-3/2 BC model sensitively depend on the selected 
value of a, and the system exhibits four kinds of behav- 
ior: For < a < 1/4, phase transition line between 
ferromagnetic F 3 / 2 and the ferromagnetic F 1 / 2 phases 
consists of only first order phase transition points and 



also it is found that a = 1/4 is independent of the co- 
ordination number. Two successive discontinuous phase 
transition lines between -F3/2 and F± phases and between 
Fx and -F1/2 phases occur for the interval 1/4 < a < 1. 
For a = 1, a discontinuous phase transition line exists 
separating F3/2 and F\ phases. On the other hand, two 
discontinuous phase transition lines emerge for a > 1, 
the aforementioned transitions are between F\ and F3/2 
phases, and between F3/2 and F\ phases. 



A. Magnetic Gibbs energy production, kinetic 
equations and relaxation time 

In order to understand the relaxation behavior of the 
spin-3/2 BC model with quenched random crystal field, 
one supposes a small uniform external magnetic field is 
applied along z axis. When a system is disturbed from 
its equilibrium state, quantities decay rapidly to their 
equilibrium values [45|. In that case, new energy form of 
the considered magnetic system in the neighborhood of 
equilibrium state can be expressed in the following: 



g(m, t, d, b) = g°(m Ql t, d, b ) + Ag, 



(6) 



where g°(mo, t, d, bo) is equilibrium magnetic Gibbs en- 
ergy in the absence of 6, and Ag is the energy variation 
of the system and it is expressed as follows: 



Ag = iAi(m-m ) 2 +A 2 (6-6 ) 

+ iA 3 (6 - b ) 2 + A 4 (m - m )(b - b ), 



(7) 



where the coefficients Ai(i — l,.-,4) can be easily cal- 
culated by using Taylor expansion and are given in Ap- 
pendix. According to theory of irreversible thermody- 
namics, the generalized force (X m ), which conjugates the 
current (m) , is derived by differentiating Ag with respect 
to m — toq : 



dAg 



d(m — mo) 



A 1 (m-m )+A 4 (6-& )- (8) 



The magnetic system will be relaxed back to its equilib- 
rium state with a small deviation in the external mag- 
netic field strength, and related process can be investi- 
gated by using linear response theory because the devi- 
ation of the applied external magnetic field is so small. 
The linear relation between the generalized flux and force 
can be written in terms of phenomenological coefficient 
as: 

m = jX m , (9) 



m 



7 Ai(m-m ) ± 7 A 4 (6-6 ), (10) 



where 7 is a kinetic rate coefficient. By assuming a simple 
form for the solution m — mo ~ exp(— t/r), the linearized 
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equation of motion has been solved and the relaxation 
time can be found in the following form: 



By solving Eq. (fTT]) numerically for a combination of 
Hamiltonian parameters we obtained the dynamic evo- 
lution of the relaxation time of the spin-3/2 BC model 
with quenched random crystal field. 



B. Derivation of kinetic equations leading to 
complex magnetic susceptibility 

When a ferromagnetic material is subject to a period- 
ically varying time-dependent magnetic held at an an- 
gular frequency w, all quantities will oscillate near the 
equilibrium state at this same angular frequency: 



III. RESULTS AND DISCUSSION 

In this section, we will focus our attention on the 
nonequilibrium dynamics properties of the spin-3/2 BC 
model with quenched crystal held. This section is divided 
into two parts as follows: In section ITlIAl we discuss how 
the quenched random crystal field affect the relaxation 
time behavior in the vicinity of order-disorder transition 
temperature as well as near the isolated critical point. In 
addition, we examine the crystal field variation of the re- 
laxation time for selected Hamiltonian parameters, and 
also we calculate dynamical mean field critical exponents 
in order to formulate the critical behavior of relaxation 
time of the considered system. Real and imaginary parts 
of the complex magnetic susceptibility are analyzed near 
the second-order phase transition temperature and also 
near the isolated critical point in section UlI Bl As a final 
investigation, dynamical mean field critical exponents are 
calculated for magnetic dispersion and absorption factors 
for both low and high frequency regions. 



m — mo = mi exp (iuii), b — bo = b\ exp (iujt). (12) 

One can easily obtain the following form substituting Eq. 
(fl"2j) into the kinetic equation given by Eq. (fT0|) : 



(13) 



mi (7A1 — iuj) + 7A4&1 = 0. 

Solving Eq. (fT3|) for m\jb\ gives: 

mi 7A4 
b\ iuj — 7A1 



(14) 



Eq. (Ti"4l can be used to calculate the complex mag- 
netic susceptibility of the considered system x(w) and the 
spin-3/2 BC model induced magnetization (total induced 
magnetic moment per unit volume) is given by: 



m — moo = Be (mi exp(zwi)) 



(15) 



where moo is the magnetization induced by a time depen- 
dent oscillating magnetic field. Also, by benefiting from 
definition, the expression for x( w ) may be written: 

m — nioo = Re (x(w)&i exp(iwt)) , (16) 

where x(w) = \ ( w ) — *X i s the complex magnetic 
susceptibility and its real and imaginary parts are called 
as magnetic dispersion and absorption factors, respec- 
tively. Finally the magnetic dispersion and absorption 
factors can be found in the following forms: 



x M 



7 A 4T 
1 + W 2 T 2 ' 

7A4T 2 W 
1 + UJ 2 T 2 ' 



(17) 



A. Relaxation time treatment of the spin-3/2 BC 
model with quenched crystal field 

Fig. QJa) represents the phase diagram in the reduced 
crystal field and reduced temperature (d, t) plane for 
a = 0.05. We should note that this type of phase dia- 
gram can be observed by selecting any value in the range 
of < a < 1/4. Here, the full and dotted lines denote the 
continuous and discontinuous phase transition points, re- 
spectively, and also the solid circle indicates the isolated 
critical point. One can clearly see from the Fig. Ufa) 
that the ferromagnetic phase separates the paramagnetic 
phase with a phase boundary line containing continuous 
phase transition points, and also phase transition line 
between ferromagnetic F 3 / 2 and the ferromagnetic Fx/2 
phases consists of only first order phase transition points 
at the relatively lower temperatures. Fig. [IJb) shows the 
thermal variation of the relaxation time r in the neigh- 
borhood of order-disorder phase transition temperatures 
for selected values of d, 7 and a. In this figure, the 
solid and dashed dotted curves correspond to values of 
d = —3.0 and —4.0 for a = 0.05 and 7 = —0.1, respec- 
tively. The vertical dashed lines refer to the phase tran- 
sition temperatures for each value of the reduced crystal 
field. As seen in Fig. [IJb) r grows rapidly with increasing 
temperature and diverges as the temperature approaches 
the order-disorder phase-transition points for both crys- 
tal field values d = —4.0 and —3.0. In order to under- 
stand the influences of the Onsager rate coefficient on 
the relaxation time, we give the Fig. HJc) for varying 
7 values, and selected values of a = 0.05 and d = 4.0. 
On both sides of the critical temperature, illustrated by 
a vertical dashed line, the relaxation time is divergent 
for all 7 values. It is worthwhile to stress that even 
though the divergence of r gets pushed away from the 
critical temperature when the absolute value of the On- 
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sager rate coefficient decreases, there exists a qualitative 
consistency between their thermal variations. Aforemen- 
tioned behaviors are in accordance with the previously 
published studies (l9M2^. On the other hand, Fig. [Tfd) 
corresponds to crystal field dependence of r near the first 
order phase transition point between ferromagnetic -F3/2 
and F1/2 phases for values of a = 0.05, t = 0.6 and 
7 = — 0.1. In contrary to treatment of the second order 
phase transition, r exhibits a finite discontinuity at the 
first order phase transition point. 

In the following analysis, let us focus on the influences 
of a parameter on the relaxation dynamics of the consid- 
ered magnetic system. In Fig. |5Ja), we plot the phase 
diagram in (d, t) plane with selected value of a — 0.3. It 
can be easily seen from the Fig. [U[a) that there exists 
two successive first order phase transitions at the rela- 
tively low temperature. The first transition is from fer- 
romagnetic -F3/2 phase to ferromagnetic F% phase while 
the second transition is between ferromagnetic F\ and 
F1/2 phases. As mentioned in section [TT1 similar type of 
phase diagrams can be obtained by focusing the range of 
1/4 < a < 1. On the other side, Fig. HJb) demonstrates 
that the r exhibits two successive infinite discontinuous 
for values of a — 0.3, t — 0.47 and 7 = —0.1 in which 
the first one is between ferromagnetic F 3 / 2 and ferromag- 
netic F\ phases and second one is between ferromagnetic 
Fi and ferromagnetic Fi/2- 

In Fig. 03a), for a bigger a parameter, for example 
a = 1.0, a discontinuous phase transition line exists sep- 
arating F3/2 and F\ phases. The crystal field variation of 
the r is shown in Fig. EHb) corresponding to the same a 
parameter for values of t = 0.4 and 7 = —0.1. It can be 
clearly seen from the Fig. EHb) that as the crystal field in- 
creases, r makes a sharp cusp at the critical crystal field 
value. As another characteristic feature of the considered 
system we have examined the crystal field variation of r 
in the neighborhood of isolated critical point. As it was 
expected, in Fig. EJc), we found that r shows a diver- 
gence treatment for selected values of parameters such as 
a = 1.0, t = 0.51 and 7 = -0.1. 

The divergence of the relaxation time on both sides 
of the critical point can be described by the critical ex- 
ponents. Thus, we may assume that, for temperatures 
lower than the transition temperatures, r obeys the law 
of the form r(i) ss (t c — t) x and values of the A can be 
obtained from the plot of log(r) — log(e), where e = 
[45| . In Fig. EJa), we have plotted r versus dimensionless 
parameter (e) in the neighborhood of continuous phase 
transition point, and we have obtained only one linear 
part on this plot for r. It is found that the slope of 
the curve equals to — 1.004 ± 0.001 indicating that r di- 
verges with a mean field exponent of —1.004 ±0.001 as 
(t c — t) _1 004±0 001 for the selected values of parameters. 
In addition to these, a similar analyze has been done for 
isolated critical point in Fig. Efib). It is deduced from the 
log-log plot of r versus dimensionless parameter that the 
slope of the curve equals to — 0.997±0.001 which indicates 
that the r of the spin-3/2 BC model with quenched ran- 



dom crystal field diverges with a mean-field exponent of 
-0.997±0.001 as (t c -t)~°- 997±omi for values of d = 1.25, 
a = 1.0 and 7 = —1.0. 

It can be seen clearly from the Fig. [5ja) that the mag- 
netic system exhibits reentrant behavior at the relatively 
higher temperature regions for a — 2.5. This behavior 
originates from the competition between negative and 
positive values of crystal field interactions. And also, 
two discontinuous phase transition lines occur at the low 
temperature regions. The first one is between the fer- 
romagnetic F\ and ferromagnetic F 3 / 2 phases, and the 
second one is from the ferromagnetic F 3 / 2 to F\ phases. 
In order to show the reentrant behavior on the relax- 
ation times, we give the Fig. [SJb) for varying values of 
7 and selected values of t, and d. As was expected, two 
successive divergences indicating two successive second 
order phase transition points exist in Fig. [5jb) for values 
of a = 2.5, t — 4.5 where the first one is from ordered 
phase to disordered phase while the other is between the 
disordered and ordered phases, respectively. It is nec- 
essary to say that as the absolute value of the Onsager 
rate coefficient decreases (7), divergence of r gets pushed 
away from the critical temperature. As a final investiga- 
tion, we focused on the low temperature region and we 
gave the crystal field variation of r in Fig. EJc). One can 
see from the figure Fig. [5jd), r shows two infinite dis- 
continuous corresponding to phase transition points, as 
crystal field interaction increases for selected values such 
as a — 2.5, t = 0.4 and 7 = -0.1. 



B. Complex magnetic susceptibility treatment of 
the spin-3/2 BC model with quenched crystal field 

Figs, inja-d) illustrate the real and imaginary parts 
of the ac susceptibility in the vicinity of isolated criti- 
cal point for both low and high frequency regions and 
values of d — 1.25, a = 1.0 and 7 = —0.1. In Figs. [HJa- 
b) the temperature variations of the magnetic dispersion 
factors are given, and, one can easily see from the figures 
that the x increases rapidly with increasing tempera- 
ture and exhibits a divergence near the isolated critical 
point for the low frequency values. On the other hand, 
the x shows a shallow dip at the isolated critical point 
for the relatively higher applied field frequencies. In ad- 
dition, two frequency dependent local maxima occur at 
before and after the isolated critical point, respectively 
(see Fig. EJb)). From the Figs. EJc-d) we see that the 
X exhibits a divergence at the isolated critical point for 
the low frequency regions, while it displays a local max- 
ima for the relatively higher applied field frequencies. By 
comparing figures EJa) and (c) one can see that the di- 
vergence treatment of \ is independent of the applied 
field frequency whereas divergence of \ sensitively de- 
pends on the frequency. On the other hand, it is possible 
to say that the similar ac magnetic susceptibility behav- 
ior can be observed by focusing any second order phase 
transition point. 
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In the following analysis, let us investigate the reen- 
trant behavior properties on the magnetic dispersion and 
absorption factors corresponding to the phase diagram 
in Fig. El[a). In Fig. (7{a) we represent the crystal field 
variation of the magnetic absorption factor for selected 
values of t, 7 and a with different applied field frequen- 
cies in the low frequency region. One can observe from 
the figure that the magnetic system exhibits two succes- 
sive divergences behavior indicating two successive sec- 
ond order phase transition points where the first one is 
from ordered phase to disordered phase while the other 
is between the disordered and ordered phases, respec- 
tively. On the other hand, as seen from figure [T^b), the 
magnetic absorption factor has two frequency dependent 
shallow dips at the second order phase transition points 
in the high frequency region. In addition, it is useful to 
mention that as we increase the applied field frequency, 
the maximum located at a finite crystal field value in the 
ordered and disordered phases decreases. 

As a final investigation, dynamical mean field criti- 
cal exponents for isolated critical point are calculated 
for magnetic dispersion and absorption factors for both 
low and high frequency regions in Fig. [5] It is known 
that the critical behaviors of x' an d x" exponents in the 
neighborhood of critical point as well as isolated critical 
point are characterized by the critical exponents. Thus, 
we may assume that, for example, when t — > ti c , y' and 
x" follow scaling laws of the form x' ^ |t — t; c | A and 
x" ~ |t — tj C I A , respectively. Here, A', A" and t; c are dy- 
namical mean field critical exponents of the magnetic dis- 
persion and absorption factors and isolated critical point 
of the considered system, respectively. In order to cal- 
culate the value of the dynamical critical exponents one 
should sketch double logarithmic plot of the related quan- 
tity versus log(I — t/t; c ) and find the slope of the linear 
part of this curve. This definition is valid for all values 
of the exponents where the negative values correspond 
to the divergence of the related quantity as 1 — t/ti c goes 
to zero and the positive values correspond to logarith- 
mic divergence, cusps or jump singularities [HI, [45[ . Fig. 
IHJa) represents the the plots of log(x') and log(x") ver- 
sus log(l — t/ti c ) at several values of applied field fre- 
quency corresponding to high region in the vicinity of 
isolated critical point for values of a = f.O, d = 1.25 
and 7 = —0.1. We should mention here that we have 
found only one linear part on each of these log — log plots, 
namely, one value for A', and one value for A", respec- 
tively. It is clearly seen from the figure that the behavior 
of log(x') and log(x") are dependent the applied field 
frequency, and the slopes of the log(x') and log(x") are 
equal to 1.0 and 0, respectively. Namely, although the 
behaviors of the magnetic dispersion and absorption fac- 
tors depend on the field frequency in the high frequency 
region, their slopes are not dependent the uj. This scaling 
form (A' = 1.0) verifies the convergence of the x' to zero 
in the high-frequency region in the neighborhood of iso- 
lated critical point (see Fig. |UJb)). In addition, scaling 
form (A" = 0.0) of the x" implies the existence a local 



maxima in the high frequency region, in other words, it is 
in accordance with Fig. H](d). On the other hand, in Fig. 
[5](b), we depict the the plots of log(x') and log(x") versus 
log(l — t/ti C ) at several values of applied field frequency 
corresponding to low region in the vicinity of isolated 
critical point for the same values with Fig. [5] (a). One 
can figure out from Fig[S] (b) that the logarithmic be- 
havior of magnetic dispersion factor does not depend on 
the applied field frequency whereas the magnetic absorp- 
tion factor changes with ui in the low frequency region. 
The mean field dynamic critical exponent of the log(x') 
and log(x") are equal to —1 and —2, respectively. It is 
obvious that these scaling forms (A' = —1.0, A" = —2.0) 
correspond to divergence behavior in the vicinity of iso- 
lated critical point in the low frequency region and for 
selected system parameters. 



IV. CONCLUDING REMARKS 

In this study, we have studied the nonequilibrium dy- 
namics of the spin-3/2 BC model with quenched random 
crystal field by means of Onsager irreversible thermo- 
dynamics. In order to obtain the Gibbs free energy of 
the considered disordered system, we used the MFA, and 
by benefiting from the minimization condition of Gibbs 
free energy, we found the mean field equation of state 
and we represented the phase diagram in the (d, t) plane 
for different a parameters. After that, we have calcu- 
lated the Gibbs free-energy production in the irreversible 
process and by means of this generalized force and flux 
are defined in irreversible thermodynamics limit. Then, 
the kinetic equation for the magnetization is obtained by 
using linear response theory. From the solution of the 
kinetic equation in the vicinity of equilibrium states, a 
relaxation time describing the nonequilibrium dynamics 
of the considered system is obtained. In order to under- 
stand the critical behavior of the relaxation phenomena, 
the temperature and also crystal field dependencies of 
the relaxation time are determined in the vicinity of the 
phase-transition points. We found that r grows rapidly 
with increasing temperature and diverges as the tem- 
perature approaches the order-disorder phase-transition 
points as well as in the vicinity of isolated critical point 
for selected Hamiltonian parameters and Onsager rate 
coefficients. And also we calculated the dynamical mean 
field critical exponents for order-disorder phase transition 
and isolated critical points, respectively and it is found 
that the slopes of the curves are equal to — 1.004 ± 0.001 
and —0.997 ± 0.001 indicating that r diverges with in- 
creasing temperature. In contrary to treatment of the 
second order phase transition, t exhibits a finite discon- 
tinuity at the first order phase transition point for fixed 
temperature value. In addition, as mentioned above we 
have touched upon the influences of Onsager rate coef- 
ficient on the r behavior, and it is observed that even 
though divergence of r gets pushed away from the criti- 
cal temperature when the absolute value of the Onsager 



rate coefficient decreases (7), there exists a qualitative 
consistency between their thermal variations. Aforemen- 
tioned behaviors are in accordance with the previously 
published studies (HHUl. 

Moreover, we have analyzed the complex magnetic 
susceptibility treatment of the spin-3/2 BC model with 
quenched random crystal field under the time dependent 
oscillating magnetic field. It is worthwhile to mention 
that the amplitude of applied magnetic field is so small 
that we made use of linear response theory in study- 
ing the magnetic relaxation process. We especially have 
touched upon the thermal variation of complex magnetic 
susceptibility near the isolated critical point for varying 
applied field frequencies. We found that the magnetic 
dispersion and absorption factors grow drastically and di- 
verges with increasing temperature in the low frequency 
regions (see Figs. HJa) and (c)). In addition, we can say 
that the divergence treatment of \ is independent of the 
field frequency whereas divergence of \ sensitively de- 
pends on the frequency. On the other hand, the x an d 
X show a shallow dip and a local maxima at the iso- 
lated critical point for the relatively higher applied field 
frequencies, respectively. It is also observed that two 
frequency dependent local maxima occur at before and 
after the isolated critical point, respectively. Moreover, 
we have investigated the crystal field dependence of the 
complex magnetic susceptibility in the neighborhood of 
reentrant region for selected Hamiltonian parameters and 
we found that the magnetic system exhibits two succes- 
sive divergences behavior indicating two successive sec- 
ond order phase transition points in the low frequency 
region. For higher frequency region, the magnetic ab- 
sorption factor has two frequency dependent shallow dips 
at the second order phase transition points. In addition, 
it is useful to mention that as we increase the applied field 
frequency, the maximum located at a finite crystal field 
value in the ordered and disordered phases decreases. As 
a final investigation we have calculated the dynamical 
mean field critical exponents of complex magnetic sus- 
ceptibility for isolated critical point. We also mention 
that a qualitatively good agreement exists between our 
findings and previously published studies [H, [ljl |4(J . 
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V. APPENDIX 



The coefficients Aj(i = 1, ..,4) can be easily calculated 
by using Taylor expansion and are given as follows: 
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FIG. 1: (a) Phase diagram of the spin-3/2 BC model with quenched random crystal field for a — 0.05. The full and dotted 
lines denote the continuous and discontinuous phase transitions, respectively. The solid circle represents the isolated critical 
point, (b) Thermal variation of the relaxation time for different values of the crystal fields near the critical point, (c) Effects 
of the 7 parameters on the relaxation time in the neighborhood of critical point, (d) Crystal field variation of the relaxation 
time in the neighborhood of discontinuous phase transition. 




FIG. 2: (a) Phase diagram of the spin-3/2 BC model with quenched random crystal field for a = 0.30. The full and dotted lines 
denote the continuous and discontinuous phase transition points respectively. The solid circles represent the isolated critical 
points, (b) Crystal field variation of the relaxation time in the neighborhood of discontinuous phase transitions. 
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FIG. 3: (a) Phase diagram of the spin-3/2 BC model with quenched random crystal field for a = 1.0. The full and dotted lines 
denote the continuous and discontinuous phase transition points respectively. The solid circle represents the isolated critical 
point. Crystal field variations of the relaxation time in the vicinity of discontinuous (b) and isolated critical point (c). 
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FIG. 4: Log-log plot of relaxation time vs. reduced parameter, and in the inset, relaxation time vs. temperature are shown 
near the continuous phase transition point (a) and near the isolated critical point (b). 
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FIG. 5: (a) Phase diagram of the spin-3/2 BC model with quenched random crystal field for a = 2.5. The full and dotted lines 
denote the continuous and discontinuous phase transition points, respectively. The solid circles represent the isolated critical 
points, (b) Effect of the kinetic rate coefficient on the reentrant behavior at t = 5.6 (c) Effect of the kinetic rate coefficient 
on the reentrant behavior at t = 4.5 (d) Crystal field variation of the relaxation time in the neighborhood discontinuous phase 
transitions. 
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FIG. 6: Temperature variations of magnetic dispersion and absorption factors for the low (a and c) and high frequency (b and 
d) region at several values of the frequency (cj) in the vicinity of isolated critical point while d = 1.25, a — 1.0 and 7 = —0.1. 
The arrow refers to the isolated critical point. 
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FIG. 7: Crystal field variations of magnetic dispersion and absorption factors for the low (a and c) and high frequency (b and 
d) regions at several values of the frequency (w) in the vicinity of reentrant region while t = 4.5, a = 2.5 and 7 = —0.1. The 
arrows refer to the critical crystal field values. 




FIG. 8: Double logarithmic (log-log) plots of magnetic dispersion and absorption factors versus (1 — t/ti c ) at several values in 
the high (a) and low (b) frequency regions while d — 1.25, a = 1.0 and 7 = —0.1. Solid and dashed lines denote the log(x') 
and log(x"); respectively. 



